Mode-coupling theory predictions for the dynamical transitions of the partly pinned 

fluid systems 
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The predictions of the mode-couphng theory (MCT) for the dynamical arrest scenarios in a partly 
pinned (PP) fluid system are reported. The corresponding dynamical phase diagram is found to be 
very similar to that of a related quenched-annealed (QA) system. The only signiflcant qualitative 
difference lies in the shape of the diffusion-localization lines at high matrix densities, with a re-entry 
phenomenon for the PP system but not for the QA model, in full agreement with recent computer 
simulation results. This flnding clearly lends support to the predictive power of the MCT for fluid- 
matrix systems. Finally, the predictions of the MCT are shown to be in stark contrast with those 
of the random first-order transition theory. The PP systems are thus confirmed as very promising 
models for tests of theories of the glass transition. 



A partly pinned (PP) fluid system is a model of a fluid 
in contact with a disordered substrate which is obtained 
by instantaneously arresting the motion of a fraction of 
the particles in an equilibrium bulk system and letting 
the remaining mobile fraction evolve under the influ- 
ence of the static random environment generated by the 
pinned particles. The pinning process can be homoge- 
neous, if, for instance, the pinned particles are randomly 
chosen in the whole volume of the bulk system, or het- 
erogeneous, if the arrested particles are those located in 
a predefined region of space. In the first case, one gets a 
model of a fluid adsorbed in a statistically homogeneous 
disordered porous solid; in the second one, a model of a 
fiuid confined by amorphous rough walls. 

Because of their peculiar preparation process from an 
equilibrium bulk fluid in which the future mobile and 
pinned components influence each other, these systems 
display a number of interesting features. For instance, 
one can show that the realization-averaged configura- 
tional properties of a PP system exactly match those of 
the bulk fluid from which it is prepared [l|, Stated 
otherwise, the pinning step does not alter the average 
configurational properties of the system. This obviously 
represents a major simplification when studying these 
quenched-disordered models, which, over the years, have 
been considered in various fields of liquid state theory. 

In the theory of adsorption in disordered porous solids, 
they appear as special cases of fluids adsorbed in depleted 
or templated matrices. Indeed, the first class corresponds 
to models of porous solids obtained by freezing equilib- 
rium configurations of a fluid and removing at random 
a fraction of its particles 0, |3j , the second one to mod- 
els obtained by freezing equilibrium configurations of a 
binary mixture and removing one component called the 
template Thus, it is clear that a homogeneous 

PP system corresponds to a depleted or templated fluid- 
matrix system in which the depleted or template par- 
ticles are reinjected as the confined fiuid 2]. Both the 
depleted and templated models are extensions of the con- 
cept of quenched-annealed (QA) system, first introduced 
by Madden and Glandt Q , in which the porous matrix is 
obtained by simply quenching equilibrium configurations 



of a fiuid without any subsequent particle removal. The 
properties of these different models have been compared 
by Van Tassel et al. j^-Q. An interesting finding that 
will be useful in the following is that, for a given matrix 
density, the depleted and templated matrices generically 
show stronger correlations and a more open accessible 
volume than a simple quenched matrix [3|-|5|. 

In computational studies of the dynamics of fiuids in 
confinement have also regularly appeared PP systems 
[H [l0l - [2H . Here, the use of these models can turn out 
highly advantageous. Indeed, the configurations of the 
fiuid particles obtained immediately after the pinning 
step are automatically equilibrium configurations by con- 
struction Hence, an acceptable initial configura- 
tion is always readily available, while finding one can be 
a very nontrivial task for other models (see Ref. [l^ for 
the case of QA systems) , and it does not have to be equi- 
librated [ll.[l3|. This is clearly very interesting for studies 
of dense and/or glass-forming systems. 

Following insightful pioneering works where PP sys- 
tems had been investigated in order to measure the spa- 
tial extent of dynamical correlations in confined glass- 
forming liquids [ll, [T2| - [l^ , it has been realized that these 
systems could also be used to probe the existence of non- 
trivial static correlations in bulk glassy liquids. The key 
here is to consider the pinning process as a constraint im- 
posed on the fluid and to measure how the conflgurations 
of the constrained particles influence the accessible states 
of the free ones, via the computation of point-to-set cor- 
relations |23l - [27| . for instance. Numerous studies of PP 
systems have recently appeared along this line [28l - [36j . 

Finally, it might be worth mentioning that an experi- 
mental nearly two-dimensional realization of such a sys- 
tem is possible, by simply squee zing a binary colloidal 
mixture between two glass plates |37l [38| . 

In this paper, we contribute to the study of the glassy 
dynamics in PP systems and report on the predictions of 
the mode-coupling theory (MCT) 39] for the dynamical 
arrest scenarios in the homogeneous case. To this end, 
we use a recent extension of the MCT to fluids imbibed in 
disordered porous solids [40l - |43j and compute the dynam- 
ical phase diagram of the pinned system deriving from 



2 



the one-component hard-sphere fluid, which, thanks to 
the overall mathematical structure of the theory and the 
robustness of the associated bifurcation schemes [sl] , can 
be expected to be a representative example. 



The motivation for this work is twofold. First, com- 
puter simulation studies of the slow dynamics of fluids 
confined in disordered porous matrices have recently ap- 
peared [19I - I22I 13 EH, allowing comparisons with the 
previous predictions of the MCT [40-43]. In many re- 
spects, they demonstrate that the MCT offers a consis- 
tent, though idealized, picture of the dynamics of these 
systems. For instance, the observed changes in the den- 
sity fluctuation relaxation pattern with increasing con- 
flnement agree well with the theoretical picture of a 
crossover from type B (bulk-like) bifurcation scenarios at 
low matrix densities to type A (Lorentz-gas-like) scenar- 
ios at high matrix densities. However, a potential issue 
with the theory is also pointed out. Indeed, while the 
MCT predicts a dynamical re-entry phenomenon at high 
matrix densities for the equisized hard-sphere QA sys- 
temic sign of it is actually visible in the simulation data 
P^I - |22| . I44 lisj. But, interestingly, Kim et al. do clearly 
flnd such a re-entry phenomenon in their results for the 
PP variant of this system (l9l-[2l| . One aim of the present 
work is to show, after some clarifications, that the MCT 
is in complete agreement with the results of the computer 
simulations with respect to this re-entry phenomenon at 
high matrix densities, even with the seemingly negative 
result obtained for the hard-sphere QA system. 



Second, the effect of pinned particles on the glass tran- 
sition has recently been investigated in the framework of 
the random first-order transition (RFOT) theory [46,. 47J 
by Cammarota and Biroli, who flnd a very interesting 
and nontrivial scenario [34j|. It is widely believed that 
the RFOT theory and the MCT have strong connec- 
tions, because they have some mathematical structures 
in common. However, actual calculations sometimes re- 
veal major discrepancies, like when the spatial dimension 
is changed |48l - [5(j | . It thus seems very natural to wonder 
how the two theories compare in the present context. As 
we will show, it turns out that the scenario predicted by 
the MCT differs very signiflcantly from the one obtained 
by Cammarota and Biroli. 



The application of the MCT scheme to fluid-matrix 
systems results in self-consistent equations for the time 
evolution of (j)q (t) , the normalized connected autocorrela- 
tion function of the fluid density fluctuations, and (t>q{t), 
the autocorrelation function of the tagged-particle den- 
sity fluctuations, at wave vector modulus q [i^ H^. In 
the infinite time limit, from which the state of the sys- 
tem can be determined, they reduce to equations for 
the nonergodicity parameters fq = limt_j.oo (f)q (t) and 



fq = limt^oo </'g(i), which read 



1-/, 
_Jl_ 



(2^ 



Vq^.k/fe/|q-k| + v;i:^A , (1) 



(1) 



(2) 



q q ■ (q - k) ,e 

,2 "^fc „2 "^|q-k| 



QC QC 

(3) 



q • k q • (q - k) i 



(2) 



q(q-k) 



Tlf 



1q-k| 



^fe^lc-ki, (4) 



q-k|' 



q (q k) 



'^Iq-kl' 



(5) 
(6) 



where rif is the fluid density, Sq the connected struc- 
ture factor, Sq the disconnected structure factor, hq 
the Fourier transform of the disconnected total pair cor- 
relation function, and Cq the Fourier transform of the 
connected direct correlation function deflned as UfCq 

1 - l/S'q [Sll-i? 

While Refs. [40|-l43j explicitly mention QA systems, no 
assumption about the statistics of the disordered solid 
is actually made during the derivation of the MCT. It 
follows that it can be applied to any fluid-matrix model, 
and in particular to homogeneous PP systems. Such a 
possibility has sometimes been challenged (T], In- 
deed, as already mentioned, the pinning protocol exactly 
preserves the conflgurational properties of the bulk fluid 
from which a PP system is prepared 3 , but it can 
have a strong effect on the dynamics [lCll - l21| . At 
flrst sight, this looks incompatible with the use of the 
MCT, which is known to take only structural quantities 
as input in order to predict the dynamics. The present 
worked example should clearly show that no such diffi- 
culty actually exists, the key being that, for fluid- matrix 
systems, one has to properly take into account the fact 
that the correlation functions generically split into con- 
nected and disconnected components in the presence of 
quenched disorder, due to broken symmetries at the mi- 
croscopic level [55|. In other words, even if this is not 
visible at the level of simple conflgurational properties, 
the structure does actually change in the pinned system 
(see the Supplementary Material of Ref. [3q] for an illus- 
tration). Similar care should be taken when discussing 
the thermodynamics of PP systems, since the thermo- 
dynamic susceptibilities such as the heat capacity and 
the compressibility are given by connected correlations 
in disordered systems [5^, ^] . 

In addition to the results for the equisized hard-sphere 
PP system, those for the analogous QA system, which 
already appeared in Refs. |40l - t43j . will be recalled for 
comparison. The difference between the two models is 
subtle. In both cases, the fluid and the matrix consist 
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FIG. 1. Mode-coupling dynamical phase diagrams of the equi- 
sized hard-sphere quenched-annealed (QA) and partly pinned 
(PP) fluid-matrix models, (pf and (pm denote the volume frac- 
tions occupied by the fluid and matrix particles, respectively. 
The dotted vertical line is a guide to emphasize the re-entrant 
behavior of the diffusion-localization line of the PP system. 



of hard-sphere particles of the same diameter and the 
control parameters are the volume fractions occupied by 
the fluid and matrix particles, denoted by <p{ and 4>„ij 
respectively. Only the preparation protocol is slightly 
modified. In the QA system, the configurations of the 
matrix particles are drawn from those of an equilibrated 
hard-sphere fluid of compacity (j)^ , then a fiuid of mobile 
particles with compacity (f>{ is inserted in the obtained 
disordered samples. In the PP system, a hard-sphere 
fluid of compacity (pf+ipYri is first equilibrated, then a frac- 
tion X — (j)ai/i<Pi + 0m) of the particles is pinned down 
to form the porous matrix of compacity 0ni, while the 
remaining ones become the mobile fluid component with 
compacity 0f. So, while the matrix in the QA system is 
prepared independently of the fluid that will be adsorbed 
in it, the matrix in the PP system is prepared under its 
direct influence. This is the reason for the generically 
stronger correlations and more open accessible volume of 
the matrix in the PP system [j-Q. Note that by con- 
struction, both models coincide when 0f or 0^ vanish. 

The required structural input is computed for both sys - 
tems within the Percus-Yevick approximation (5ll - [53l[57| . 
which leads to analytic expressions for the hard-sphere 
PP system, as it was recently realized [2] • The numerical 
procedure used to solve the above equations is described 
in Refs. dJE^, to which the interested reader is referred 
for technical details. 

The resulting phase diagrams are plotted in Fig. [T] 
Both have much in common. Three phases are found: 
fluid when /g = /| = is the only solution, localized 
when fq = and fq ^ 0, and glassy when fq ^ and 
fg 0. The systems enter the glassy phase by crossing 



an ideal glass transition line, which corresponds to the 
freezing of the collective dynamics and along which the 
bifurcation scenario changes from type B to type A as 
the confinement gets stronger. Furthermore, this line is 
preceded at high enough matrix densities by a diffusion- 
localization transition line associated with a continuous 
arrest of the self dynamics only. For the PP system, 
there might be a very short glass-glass transition line at 
the junction between the type A and B glass transition 
branches, but given its very limited extent and the possi- 
ble sources of inaccuracies in the numerical calculations, 
such as the discretization of the integrals in Eqs. ([TJ and 
([2]), no definite conclusion can be reached at this point. 
In practice, both systems can be said to have phase dia- 
grams with the same topology. 

Now considering the shape of the phase diagrams, an- 
other shared feature of both models is a significantly re- 
entrant ideal glass transition line, the effect being more 
pronounced for the PP system. But, and quite remark- 
ably, the MCT delivers contrasting predictions with re- 
spect to the behavior of the diffusion-localization tran- 
sition lines. Indeed, no re-entry phenomenon is visible 
in Fig. [1] for the QA system, while one is clearly present 
for the PP system For both types of transition, the 
overall trend, which is that the ergodicity-breaking events 
occur at higher densities in the PP system, is perfectly in 
line with the fact that it is characterized by a more open 
accessible volume [^-Q , thus confirming an argument by 
Kim et al. 

The diffusion-localization transition lines are precisely 
those investigated by Kim et al. in their simulations |19l - 
[2l| . Indeed, the criterion they use to define the dynam- 
ically arrested phase is based on the mean-squared dis- 
placement, a self dynamical quantity. This is confirmed 
by the careful analysis of the interplay between the col- 
lective and self dynamics led by Kurzidim et al. [1^ . 
Hence, the conclusion: the MCT does correctly capture 
the fact that the subtle structural differences between the 
QA and PP models result in different re-entry behaviors 
of their self dynamics in the high matrix density reg ime, 
in full agreement with the simulation results |19l - [2H . 

Unfortunately, the situation is less clear with respect to 
the collective dynamics, which has been studied in detail 
for the QA system only 13] • Indeed, in qualitative 
agreement with the prediction of distinct transition lines 
for the self and collective dynamics at moderate and high 
matrix densities, a wide separation of time scales between 
the two dynamics is observed in simulations, but with no 
sign of a re-entrant collective dynamics. The problem 
might lie in the theory or in the simulations. Indeed, on 
the one hand, it is possible that the standard formulation 
of the MCT based on equilibrium quantities breaks down 
if the particles cannot redistribute themselves across the 
system, as it is precisely the case in a localized state. 
On the other hand, the simulations, which are very diffi- 
cult in this density regime, do not seem consistent with 
the expectation that the collective and self dynamics co- 
incide in the limit of vanishing fiuid density. So, more 
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FIG. 2. Mode-coupling dynamical phase diagram of the equi- 
sized hard-sphere partly pinned fluid-matrix model. (j>{ -\- 
is the total volume fraction occupied by both the fluid and 
matrix particles, x — tjjm/ (<?!>f -I- 0m) is the pinning fraction. 



work is definitely needed to clarify this point. Here, the 
PP systems could be very useful, thanks to the possibil- 
ity to easily generate equilibrated samples even at high 
densities. 

We now turn to the comparison with the predictions of 
the RFOT theory [sj]. For this, it is more convenient to 
consider the dynamical phase diagram of the PP system 
in the plane defined by the total volume fraction (/)[ + (ji-^^ 
and the pinning fraction a;, as in Fig. [5J In passing, we 
note that, thanks to this representation, the total den- 
sities at the various transitions are clearly shown to be 
decreasing functions of x, in line with the rather natural 
expectation that pinning down particles slows down the 
dynamics. 

Within the RFOT theory, the PP systems are found 
to occupy a very special place. Indeed, at different levels 
of calculation (mean-field or real-space renormalization 
group analysis), the theory consistently predicts ideal 
glass transition lines starting at x = (bulk limit) and 
ending at some critical value < 1. This means that 
one can go from the fluid to the glassy phases without 
crossing any transition, simply by moving around this 
terminal point as one can do with the critical point of 
the liquid-gas transition. However, this is an isolated 
case. For other types of random environments, it is pre- 
dicted that, if the disorder becomes strong enough, the 
glass transition does not disappear, but turns continuous, 
very much like in the MCT with the type A bifurcations, 
and it is then impossible to reach the glassy phase with- 
out crossing a transition. 

So, while the predictions of the MCT and the RFOT 
theory tend to agree for generic cases of fluids in disor- 



dered environments, they are in stark contrast for the 
PP systems. Indeed, within the MCT, the PP systems 
appear as rather ordinary fluid-matrix systems, as illus- 
trated by Fig.[l] and the transition lines run up to a: = 1, 
so that the system always has to cross a dynamical tran- 
sition line in order to enter the glassy phase, as shown 
in Fig. [21 The PP systems are thus a clear example of 
diverging predictions between the two considered theo- 
retical schemes and, as such, would probably deserve 
much attention in the future. We note that, based on 
the presently available data and, in particular, on the 
joint study of QA and PP systems by Kim et al. [l9l - 
no obvious inconsistency with the MCT scenario has 
appeared yet, but, as pointed out above, some issues are 
still pending and more work is needed. 

In summary, in this paper, the predictions of the MCT 
for the slow dynamics of a homogeneous PP system have 
been worked out and compared to previous findings for a 
related QA system pOl - |43l. t o the results of recent com- 
puter simulation studies and to the scenario de- 
rived within the framework of the RFOT theory 's?] . It 
is found that the dynamical phase diagrams of the QA 
and PP systems are very similar within MCT. The only 
significant qualitative difference is in the shape of the 
diffusion-localization lines at high matrix densities, with 
a re-entry phenomenon for the PP system but not for 
the QA model. In the light of these results, the simula- 
tion data, which could first appear as a challenge to the 
theory, are actually found to represent a nice confirma- 
tion of its predictions. Finally, the comparison between 
the MCT and the RFOT theory shows that the two ap- 
proaches make predictions that are in stark contrast. 

The latter result might have serious consequences. In- 
deed, two different points of view might be adopted on 
PP systems. On the one hand, seen from the theory 
of simple liquids and adsorption phenomena, they are 
fiuid- matrix models among others They display 

interesting and possibly useful configurational properties 
[l], m, but these are merely accidental and no peculiar 
physics beyond confinement phenomena should be ex- 
pected. This is what the MCT supports. On the other 
hand, seen from the theory of disordered systems, they 
are very special physical objects, ideally and uniquely 
suited to track hidden forms of order and the related 
phase transitions in amorphous systems [23l - l36j . This is 
what the RFOT theory illustrates. Finding which point 
of view is the most appropriate for glass-forming liquids 
could have a profound influence on future directions in 
glass transition theory. 
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